Dynamics of coupled cavity arrays embedded in a non-Markovian bath 
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In this paper, the non-Markovian dynamics of a coupled A^-cavity model is studied based on 
the quantum state diffusion (QSD) approach. The time-local Diosi-Gisin-Strunz equation and the 
corresponding master equation are derived for a coupled cavity array. Non-Markovian effects are 
studied in two examples, two-cavity and three-cavity, under different boundary conditions. We have 
shown that the environment-memory can facilitate the cat-like state transfer from one cavity to 
another in the case of a strongly non-Markovian environment. We show that the non-Markovian 
QSD approach can be a valuable tool for the dynamics of a multi-partite continuous- variable (CV) 
system. 

PACS numbers: 42.50.Lc, 03.65.Yz, 05.30.-d 



I. INTRODUCTION 

The dynamics of quantum open systems is one of the 
most important research topics in quantum optics, 
tum dissipative system and quantum information 
Open system dynamics is also essential for a deeper un- 
derstanding of decoherence and implementation of quan- 
tum control While the theory of open systems in 
the Markov regime is well understood using the standard 
Lindblad master equations or the corresponding Markov 
quantum trajectories [l3l - [l6| . the recent surge of interest 
in non-Markovian phenomena requires a set of system- 
atic theoretical tools for approaching quantum dynam- 
ics in strong non-Markovian regimes. Among many re- 
cent research efforts in developing a theoretical approach 
to investigating the non-Markovian evolution, the non- 
Markovian quantum state diffusion (QSD) equation ini- 
tially proposed by Diosi, Gisin, and Strunz [13, [ill pro- 
vides a very powerful tool to deal with quantum open sys- 
tems coupled to a bosonic environment. The Diosi-Gisin- 
Strunz QSD equation was derived directly from an un- 
derlying microscopic model irrespective of environmen- 
tal memory time and coupling strength. Several physi- 
cally interesting models have been solved using the non- 
Markovian Q^SD approach ranging from one qubit [l9| . 
two-qubit 1201 . iV-qubit systems [2l[, to multi- level sys- 
tems [13, |23| |. Moreover, recent research has extended 
the theory of QSD for a bosonic bath to an open system 
coupled to a fermionic bath [13, HI]. Furthermore, as a 
computing tool in real applications, the QSD approach 
also showed its potential value in many other interest- 
ing problems including a theory of p recise measurement 
[27[ , quantum control dynamics |28| , and quantum biol- 
ogy [2^. 

The non-Markovian QSD equation has been widely 
appHed to discrete- variable systems including qubits. 



qutrits, spins etc. However, it is also important to 
consider the non-Markovian QSD approach for multiple 
continuous-variable (CV) systems since CV systems are 
known to have many unique features leading to advan- 
tages in its application in quantum information process- 
ing (QIP) such as quantum teleportation, quantum error- 
correction codes, dense coding, and quantum computing 
[sol ]. Furthermore, CV systems are an integral part of 
quantum dissipative systems, non-equilibrium statistical 
physics and theory of quantum decoherence [^ . 

The purpose of this paper is to derive the non- 
Markovian QSD equation and the corresponding master 
equation for a coupled A^-cavity model. This is a very 
good example of using the QSD technique to solve a mul- 
tiple CV quantum open system. It should be noted that 
the coupled A''-cavity model itself is an important model 
in QIP due to its relevance to the quantum communica- 
tion network schemes 
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As simple cases of the coupled 
A^-cavity model, single-cavity and coupled two-cavity 
have been widely studied and the exact master equations 
are derived by using the path integral approach [3, [^, [13] . 
We extend the two-cavity case into a multi-cavity case 
and show that the QSD technique is a very useful tool 
in dealing with the A-cavity model. Then, based on 
the non-Markoivan QSD equations and master equations, 
we also study several interesting phenomena concerning 
quantum coherence transfer assisted by the environment- 
memory. Specifically, in the two-cavity example, we show 
that the memory-effect of highly non-Markovian environ- 
ment can be used to transfer the so-called Schrodinger 
cat-like state between the two cavities. This phenomenon 
exists only in the strongly non-Markovian environment, 
so it may be called the memory- assisted state transfer. 
Furthermore, in the case of the three coupled cavities 
model, we study the cat-like state transfer and entangle- 
ment transfer under different boundary conditions, and 
show that the dynamics of the system is strongly depen- 
dent on the boundary conditions. 

This paper is organized as follows: In Sec.|lll we derive 
the exact QSD equation and the corresponding master 
equation for coupled A^-cavity model. The finite temper- 
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ature case is subsequently discussed in Sec. IIIIl Then, in 
Sec. llVi we study a special example, the coupled two- 
cavity case. We show that the cat-like state can be 
transferred from cavity 1 to cavity 2 without any direct 
couplings. Furthermore, we show that the cat-like state 
transfer can only be observed in a highly non-Markovian 
environment. Next, in Sec. |Vl the coupled three-cavity 
example is considered. We show that the quantum dy- 
namics (cat-like state transfer or entanglement transfer) 
will be significantly modified by the boundary conditions. 
Finally, we conclude in Sec. |Vll 



II. COUPLED TV-CAVITY MODEL: 
ZERO-TEMPERATURE BATH 

We consider a coupled cavity array interacting with 
a common bosonic bath. The total Hamiltonian of the 
system plus environment can be written as, 

^^tot ^ Hs + Hb + Hint, (1) 

where the Hamiltonian of the system is described by 



N 



N 



1=1 i=l 

where uJi {i = 1,2,...N) are the frequencies for N cavi- 
ties, and Xi are the coupling constants between ith cavity 
and (i + l)th cavity. We consider two types of boundary 
conditions for this model. The first one is the periodical 
boundary condition (PBC), which means that cavities 
form a closed loop (i.e., flAr+i = ai). The other one is 
the open boundary condition (OBC). Namely, the cavity 
chain has two open ends. We assume that all the cavities 
are coupled to a common bath, which is composed of a 
set of harmonic oscillators 6 , , b^. 



(3) 



The interaction Hamiltonian between the system and the 
environment is 



(4) 



For simplicity, we assume that each kij — gjU is the 
product of two coupling constants gj and U , then 



(5) 



where the Lindblad operator L — X^iLi ho,i- In the fol- 
lowing derivation, we set all the coefficients = 1 for 
simplicity. 

The formal Diosi-Gisin-Strunz equation can be written 
as [13, 

—ipt = -iHsi^t + Lzl^t - L'^ j dsa{t,s)j^, (6) 



where a{t,s) = X^jlSjP^ bath correla- 

tion function at zero temperature (T = 0), and = 



is a complex Gaussian process satisfy- 



ing M[zt\ = M[ztZs] = 0, M[zlZs\ = a{t,s). Here 
M[-] denotes the statistical average over the classical 
noise z^. The reduced density operator of the system 
may be recovered by the statistical average from the 
quantum trajectories ipt{z*): pt = M[\ipt{z*)) {ipt(z*)\] = 

In order to establish a time-local non-Markovian QSD 
equation, we have to replace the functional derivative 
by a time-local operator 0{t,s,z*), 

SMz* 



Szt 



Oit,s,z*)Mz*)- 



(7) 



We will use the shorthand notation: O = 0{t, s, z*) if no 
confusion arises. According to the consistency condition 
^^Vi(-2*) = ^ ^V't (-2*) J the operator O can be shown 
to satisfy the following dynamic equation, 

-0=[-tH, + Lz;-L^O,0]-L^ — 0, (8) 

where O is defined as 0{t,s,z*) = 

The initial condition is 



dsa{t, s)0{t, s, z*)ds. 
given by 
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Oit,t,z*)^L^J2' 



(9) 



From Eq. ([5]), the exact O operator of this model can 
be derived as. 



N 



(10) 



Note that the O operator for this coupled cavity array 
does not contain the noise zj* . Substituting Eq. (ITU|) into 
Eq. ([U, we can show that the coefficient functions fi{t, s) 
satisfy the following equations, 

d 

glMt, s) = lUjJi{t, s) + l[Xifi+l{t, s) + Ai_i/i_i(t, s)] 



N 



[ID 



where Fi{t) = J* dsa{t, s)f,{t, s) (i = 1,2,..,N), with 
the initial conditions 



f^{t = s,s) = 1. 



(12) 



With the exact O operator given by Eq. (|T0)) . we can 
derive the exact master equation for the system from the 
exact QSD equation [l9| . 



d 



p = -t[Hs, p] + [L, pOHt)] + [d{t)p, it] 



N 



-i[Hs, P] + {J2 FI(.t)ia^,pa]] + h.c.}. (13) 
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This exact master equation is derived in a very sim- 
ple manner due to the fact that O operator is noise- 
independent. However, when the initial state of the bath 
is a finite-temperature thermal state, the operator O will 
not be noise free, and the derivation of the exact master 
equation becomes more involved as to be shown in the 
next section. 



III. COUPLED AT-CAVITY MODEL: 
FINITE-TEMPERATURE BATH 



In the case of the finite-temperature bath, we can in- 
troduce a fictitious bath and transform the finite tem- 
perature case into an effective zero temperature model 
jisl l32j | . The non-Markovian QSD equation for the finite- 
temperature case is given by 

^Vt = [-iH, + Lzt + L^wl - L^O, - L02)A, (14) 



where Hg is given in Eq. ([2]), and — /' 



/i*u'*e~"^'* are two statistically indepen- 
dent Gaussian noises with fi — ^/fii + Igi and hi = 
y/nlgi as the effective coupling constants. Here ni = 
exp(hiy }kBT)~i ' "^^^ Gaussian noises satisfy the fol- 
lowing relations, 

M[zt] = M[ztz,] = 0, M[ztz,] = a^{t,s), (15) 
M[wt] = M[wtWs\ = 0, M[wlws\ - a2{t, s). (16) 

where ai(i,s) = T.^\f^?e~"''^'"'^ and a2{t,s) = 
X]i l/iipe*''''*^*^ are correlation functions for the two ef- 
fective baths. For the finite-temperature case, the two O 
operators are defined as 



Oi{t,s,z* ,w*)i)t 

02{t,S,Z* ,W*)^t 



diptjz* ,w*) 



Sw* 



(17) 
(18) 



satisfying the following equations 



d 



— Oi = [-tH, + Lz; + L^w; - L^Oi - LO2, Oi] 
at 



— O2 = i-iH, + Lz; + L^wl - L^Oi - LO2, O2 
at 

-L^l-0,-L^02, 
owz ow* 



(20) 



where Oi — ai{t, s)Oi{t, s, z* ,w*)ds (i ~ 1,2), Ac- 
cording to Eq. ()19ll20p . the exact O operators for the 
A^-cavity model can be determined as 



N 



Oi{t,s,w*)^y pi{t,s)ai+ / q{t,s,s)wl,ds , (21) 



1=1 



02it,s,z*) = y2^^it^M+ f y{t,s,s')zl,ds', (22) 
while the coefficients satisfying the following equations 

= iiOiPiit, s) + i[Xip.i+i{t, s) + Ai_ipi_i(i, s)] 

N N 



dt 



Xi{t,s) 

■- -iuJiXi{t, s) - i[\iXi+i{t, s) + \i^iXi^i{t, s)] 
N N 

^j{t, s)P,{t) - Y ^At^ s)X,{t) - Q{t, s), (24) 



where 











N 


(25) 




N 

i=i 


(26) 


m) 


= / ai{t,s)pi{t,s)ds, 
Jo 


(27) 


Qit,s') 


= / ai{t, s)q{t, s, s)ds, 
Jo 


(28) 


X.,{t) 


= / a2{t,s)xi{t,s)ds, 
Jo 


(29) 


Y{t,s') 


= / a2{t,s)y{t,s,s')ds, 
Jo 


(30) 



with the initial conditions 

Pt{t,t) = 1, 

Xi{t,t) = 1, 

q{t,t,s')=0, 
y{t,t,s')=0, 



N 



q{t,s,t) = -^pj(i,s), 

N 

yit,s,t) = ^a;j(t,s). 



(31) 
(32) 

(33) 
(34) 

(35) 
(36) 



Since the O operators in Eq. ^111221) are noise-dependent, 
the exact master equation cannot be obtained directly 
from Eq. (jl3p . We have to derive the new master equa- 
tion in the finite temperature case. Using the so-called 
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Heisenberg approach [2^, the exact master equation is 
derived as 
d 

—p = -i[Hs,p\ 

+ {CI (t) [J2a,,J2 P"]] + ^2 (t) [E °- E 

+ CI it) [E E z^"^-] + ^4 it) [E «I ' E «^-'°] 

+ h.c.}, (37) 
where the coefficients Ci{t) (i = 1, 2, 3, 4) are given by 

C.(t)= /" ai(i,s)c,(t,s)ds (i = l,2), 
Jo 

Cj{t)= [ a2{t,s)cjit,s)ds 0- = 3,4), (38) 
Jo 

and the details of the derivation will be shown in Ap- 
pendix A. 

This is the general form of the master equation for iV- 
cavity model in the finite temperature case. This model 
serves as a good example showing that the QSD approach 
is a useful tool in deriving the non-Markovian master 
equation in both zero temperature and finite temperature 
case. In the following two sections, we will discuss several 
examples with N — 2 and = 3. 

IV. TWO COUPLED CAVITIES 

In this and the next sections, we will show several in- 
teresting numerical results based on the simple examples 
N = 2 and A = 3. The non-Markovian QSD equation 
for the model consisting of one cavity was discussed in 
[l^ , here we start from the case of two coupled cavities. 
For the case with N = 2, the system Hamiltonian and 
the Lindblad operator are given by, 

Hs — ujialai + uj2a\a2 + \{a\a2 + aia\), (39) 

L = ai + 02. (40) 

According to the general QSD equation for A-cavity 
model in Sec. HIl the non-AIarkovian QSD equation for 
the coupled two-cavity model is 

^^iPt = i-iH, + Lzl-L^O)i^u (41) 

where 0(t) = Fi(t)ai -I- F2(i)a2, the coefficients are de- 
termined by Eq. (jlip . From Eq. (jl3p . the corresponding 
master equation for the coupled two-cavity model is 

= ~iujiia\aip - pa\ai) - iuj2{a\a2P - pa\a2) 

— i\ia\a2P — pa\a2 + a|aip — pat^ai) 

+ {Fiit)iaipa\ — a\aip + aipa\ — a^aip) 

-f F2{t){a2pa\ - a\a2P + a2pa\ - a\a2p) + h.c.}. 

(42) 
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FIG. 1: (Color online) Wigner functions for the two cavities 
are plotted at cjt = 0, ojt = 11, ujt = 22, and cot = 33. The 
initial state is chosen as the cat-state l/v'Z(|l) -I- | — 1)) ® |0) . 
The parameters are uji — UJ2 = uj = 1, \ = 0, J = 0.1. 



Note that our derivations of the QSD and master equa- 
tion are independent on the choice of the environment 
correlation functions. For simplicity, all the numerical 
simulations using the QSD equation are based on the 
Ornstein-Ulenbeck noise, which is defined by ait, s) = 
2g-7|t-s|^ The Ornstein-Ulenbeck process recovers the 
Markov limit when 7 — >■ 00, ait, s) — )■ Sit, s). This prop- 
erty will allow us to observe the crossover properties from 
non-Markovian to Markov regimes. A small 7 typically 
represents a non-Markovian noise with the finite memory 
time. 

For two coupled cavities without interacting with an 
environment, one can easily see that the phenomenon 
of cat-like state transfer due to the mutual coupling be- 
tween the two cavities. In fact, if we choose the initial 
state of the first cavity as the Schrodinger's cat-like state 
IV'cat) = 1/v^(|ck) + \ ~ ct)) (-Z^ is a normalization fac- 
tor), while the second cavity is in the vacuum state, i.e., 
l^PiO)) = l/VZi\a) + \-a))i'^\0)2, the cat-like state will 
be transferred to the second cavity from the first cavity 
after time t^: |-0(i = t^)} = \0)i (81 1/Vzi\a) + \ - a))2. 
This result can be simply proved by solving the corre- 
sponding Fokker-Planck equation [33] . 

Using the exact QSD equation and master equation, 
we study several interesting properties of non-Markovian 
dynamics of the multiple cavity model. We find that 
the non-Markovian common bath can induce a mutual 
correlation between two cavities indirectly, and make the 
cat-like state hop from one cavity to another even the two 
cavities are not interacting to each other. Notably, this 
phenomenon of the memory-assisted transfer can be only 
observed in a strongly non-Markovian environment. It is 
easy to show that the Markov environment will quickly 
suppress the possibility of the state transfer due to the 
irreversible energy dissipation to the environment. In 
Fig. [U we plot the Wigner functions for the two cavi- 
ties at different time points. The regions with negative 
values reflect the quantum interference due to the co- 
herent superposition of two wave packets. In a strongly 
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non-Markovian case (7 = 0.1), the cat-like state can be 
transferred to the second cavity from the first cavity with- 
out any direct couphng between two cavities (i.e., A = 0). 
As a comparison, it is noted that when the environment 
is in a Markov regime (7 is large) , the high fidelity trans- 
fer cannot happen (the Markov case is not plotted in 
the figure). Since the dissipative process is dominant 
in the Markov limit, the decoherence process suppresses 
the transfer process swiftly. Therefore, the cat-like state 
may be destroyed by dissipation long before it can hop 
to the second cavity. The numerical results show how 
a non-Markovian environment affects the quantum state 
evolution. 
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V. THREE COUPLED CAVITIES 

A more complicated example is the case of three cou- 
pled cavities. For iV = 3, the system Hamiltonian can be 
written as 

Hs = ujia\ai + uj2a\a2 + ^30303 + Ai(a|a2 -|- aia\) 
-f \2{a\az -I- a2a\) -|- A3(a|a3 + 0103), (43) 

and the Lindblad operator is given by, 

L = ai+ 02 + 03. (44) 

This model allows us to consider two types of boundary 
conditions. In the case of the open boundary condition 
(OBC), we have Ai = A2 = A, A3 = 0. In contrast, in 
the case of the periodic boundary condition (PBC), we 
have following constraints: Ai = A2 = A3 = A. Accord- 
ing to the QSD equation for the general A^-cavity model 
discussed in Sec. [Ill the exact QSD equation is given by 



(45) 



where O — Fi{t)ai + F2{t)a2 + ^3(^)03, and the coef- 
ficients are determined by Eq. Consequently, the 
corresponding exact master equation is given by, 

• /' t t ^ ■ / t t \ 

—p = —iuji(a{aip — pa[ai) — iuj2[a'2a2P — ^0202) 

- iuj'i[a\a'ip - palas) 

— iXi{a\a2P — pa\a2 + a^aip — pa^a-i) 

- i\2{ala3p - palas + a\a2P - pa\a2) 

— i\'i{a\azp — /9a|a3 -I- a\aip — pa\ai) 
+ {-P^r (0(«iPa| - pa\ai + a2pa\ 

- pa\a2 + azpa\ - paXa^i) 

+ F2{t){a2pa\ - pa\a2 + aipa\ 

— pa\ai + 03^02 — pa\a'i) 

+ F^{t){a3pal — pala^ + aipal 



FIG. 2: (Color online) Time evolution of fidelities, (a) is 
for the case of PBC, and (b) is for the case of OBC. The 
red (solid), green (dashed), blue (dash-dotted) lines are the 
fidelities J-^^\ J-^'^\ and J-^^^ for cavities 1, 2, 3, respectively. 
The parameters are cJi = a;2 = = w = 1, A = 1, 7 = 0.2. 



OPEN BOUNDARY CONDITION 

(ot=0 (ot=1.1 [ot=2.2 











FIG. 3: (Color online) Wigner functions for three cavities in 
the case of OBC at different time points. The initial state is 
chosen as 1/\^Z{\1) -I- | — 1)) (8) |0) ® |0) . The parameters are 
oji = aj2 = 1^3 = <^ = 1, Ai = A2 = 1, A3 = 0, 7 = 0.2. 



pagCi -|- 02^03 — pala2) + h.c.}. 



(46) 



The above exact QSD equation and master equation will 
serve as a powerful tool in our study on the properties of 
the non-Markovian dynamics under two different bound- 
ary conditions. 
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PERIODICAL BOUNDARY CONDITION 

(ol=0 cot=1.02 tot=2.04 








FIG. 4: (Color online) Wigner functions for three cavities in 
the case of PBC at different time points. The initial state is 
chosen as + | — 1)) ^ |0) (8) |0) . The parameters are 

oji = CJ2 = aj3 = oj = 1, Ai = A2 = A3 = 1, 7 = 0.2. 



Cat-like state transfer under different boundary 
conditions 



We have already showed that the cat-like state 
\l\fzi\a) + I — a)) can be transferred from one cav- 
ity to another both analytically and numerically in the 
case of = 2. However, the cat-like state transfer 
can also be realized in the TV = 3 case. Moreover, 
the model of three coupled cavities allows us to con- 
sider the two different boundary conditions that will pro- 
foundly affect the effect of state transfer. In Fig. [21 
fidelities are plotted for three cavities under two dif- 
ferent boundary conditions. The fidelities are calcu- 
lated as T^^ = max{,}[(7^(f,)(0)|p(»)(t)|V'^S)(0))], where 

|^(f^)(6»)) = 1/Vz(|ae-*^) -f I - ae-*^)) is the rotated 
cat-like state, p^'^ (i) is the reduced density matrix for the 
cavity «, and T^'^^ is the fidelity for the cavity i. Since the 
cat state may rotate in time evolution no matter with or 
without an environment, so all the rotated cat-like state 
IV'cft (^)) fo'" given d can be considered as a trans- 
fered cat-like state. Therefore, we calculate the maxi- 



mum fidelity between and IV'cat (^)) possible 

B. From the curves of fidelities in Fig. [21 it is shown that 
in the case of PBC, the fidelity will increase and reach to 
1 only in the first cavity. In the case of OBC, the fidelity 
of the third cavity can be close to 1 at « 2.2, which 
means that the cat-like state is transferred into the third 
cavity. It is worthwhile to point out that, in the case 
of OBC, the cat-like state will only be transferred from 
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FIG. 5: (Color online) Time evolution of the negativity for 
all the three modes. The red (solid), green (dashed) and blue 
(short-dashed) lines are entanglement (negativity) between 
cavities 1-2, 1-3 and 2-3, respectively. The initial state is 
chosen as |V'(0)) = ^(|100) -f |010)). The parameters are 
a;i = aj2 = ix^a = cj = 1, A = 1, 7 = 0.2. 



cavity 1 to cavity 3. It means that the state of the cav- 
ity 2 will never be a cat state. Intuitively, in the case of 
OBC, two cavities at the end points of a chain are sym- 
metric while the middle cavity is in a unique position. 
This is similar to the classical model of a mechanical os- 
cillator in a string, the superposition of the oscillation 
forms a standing wave due to the boundary conditions, 
so that the middle point on the string is just a "node" , 
whose amplitude is always zero. In contrast, for the case 
of PBC, the cat-like state never transfer into other cav- 
ities since the other two cavities are in the symmetric 
positions. The state has an equal probability to be one 
of the two cavities, so the states of two cavity cannot 
be the cat state at the same time. Therefore, one will 
expect that the cat state will experience a decay-revival 
process in the initial cavity. To examine this more explic- 
itly, in Fig. [3] and Fig. [H we plot the Wigner functions 
for each cavity in the three-cavity model with OBC and 
PBC. The time evolution of the Wigner functions shows 
that the cat-like state is able to transfer into cavity 3 at 
ijjt = 2.2 in OBC case. However, in the case of PBC, one 
only observe the decay and revival of the cat-like state in 
the first cavity (initial cavity) at ujt — 2.04. Of course, 
the cat state will disappear completely in a long-time 
limit due to the dissipation effect. 



Entanglement transfer under different boundary 
conditions 



As another interesting phenomenon, let us consider the 
entanglement transfer in the cavity array. We choose 
the initial state of the three-cavity chain as |'0(O)) = 
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^(|100) + |010)) (in Fock basis), where cavity 1 and 
2 are entangled, but the combined cavity 1 and 2 are 
separable from cavity 3. In the OBC case, the initial en- 
tanglement between cavity 1 and 2 will transfer to the 
combined system of cavity 2 and 3. Specifically, in Fig. [5] 
(OBC case), the entanglement between cavity 2 and 3 
will reach the maximum value at ujt « 2.2. At the same 
time point, the red line and green line in Fig. [5] show 
that cavity 1 become separable from both cavity 2 and 
cavity 3. As a comparison, we see that, in the PEC case, 
the entanglement between cavity 2 and 3 never reach 1. 
However, PBC case is shown to be more helpful to the 
entanglement preservation. For example, the loss of en- 
tanglement degree in cavity 1 and 2 at w 4 is less than 
that in the case of OBC. 



VI. CONCLUSION 

In this paper, we have derived the time-local non- 
Markovian QSD for the coupled A^-cavity model embed- 
ded in a common bath without any approximation. Us- 
ing the non-Markovian quantum trajectory approach, we 
have successfully derived the exact master equation for 
the N coupled cavity model in both the zero and finite- 
temperature bath cases. Based on these exact evolu- 
tion equations, we studied several interesting phenomena 
in the dynamic evolution of the two-cavity and three- 
cavity cases. We studied the cat-like state transfer in 
the coupled cavity array, and showed that a highly non- 
Markovian noise can enhance the cat-like state transfer 
(called the memory-assisted cat state transfer). In the 
three-cavity example, we also studied the effect of bound- 
ary conditions on the cat-like state transfer and the en- 
tanglement transfer. These examples are used to show 
that the non-Markovian QSD approach is a highly valu- 
able tool in characterizing the dynamics of a CV system. 
It is important to note that the approach can be directly 
extended to a generic multipartite CV system with an 
arbitrary number of cavities. 
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Appendix A: Deriving the non-Markovian master 
equation in the finite temperature case 

In this section, we will provide details of deriving the 
finite temperature master equations from the exact QSD 
equations. We only consider the PBC case, and notice 
that the OBC case can be treated in a similar way. For 
simplicity, we assume that all the cavities are identical 



(Le., (£j = w, Ai = A). Following the Heisenberg approach 
|26l |32[, one can deal with the derivation of the exact 
master equation in the general case. 

In the finite temperature case [l^ 113 , the total Hamil- 
tonian in the interaction picture is given by 

i 

+ J2 Ke-'-'^'L^el + h.e^'^'Le,), (Al) 

i 

where the Hamiltonian of the system (coupled A'^-cavity) 
can be written as 

N N 

Hs^u^ alai + A ^(a|ai+i -I- aJ+^Oi). (A2) 

1=1 1=1 

In the Heisenberg picture, the evolution for any operator 
is given by, for example, 

^d,(t) = ~iUr'[d^, H[ll{t)]Ut = -ihe'^^'L (t) . (A3) 

Integrating both sides of the equation and noticing the 
initial condition di{0) = d^, we obtain 

d,{t) =d,- if, f L{s)e"'^'ds, (A4) 

similarly, we also have 

e,{t) = d,-ih* I L\s)e-'''^'ds, (A5) 

dUs) = dlit) - if* f Lt(s')e— '^'ds', (A6) 

J S 

e\{s) = e\{t) - ih, J* L{s')e'''^'' ds' , (A7) 

For a single trajectory, the time evolution of the state 
of the system is governed by the stochastic propagator of 
the system G{t,z*,w*) = {z\{w\Ut\0) as 

\i;{t,z*,w*))^G{t,z*,w*)\iJsm)- (A8) 

Applying Schrodingcr equation to the whole system and 
noticing Eq. (|A4IIA7p . we wiU find G{t,z*,w*) satisfies 
an equation similar to the QSD equation (I14p . 

^^G{t,z*,w*) 

= -iHsG + Lz;G - [ dsai{t,s){z\{w\UtL{s)\0) 
Jo 

+ L^wlG~Lf dsa2{t,s){z\{w\UtL''{s)\0}, (A9) 
Jo 
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where ai{t,s) = Ei l/^Pe-*''^*"'^ "2(^,5) = 
Ei From this dynamic equation, we can 

also define the O operators from the stochastic propaga- 
tor 

{z\{w\UtL{s)\0) =Oi{t,s,z*,w*)G{t,z*,w*), (AlO) 

(z|(u)|t/tL^(s)|0) = 02{t,s,z*,w*)G{t,z*,w*). (All) 

Instead of deriving the non-linear differential equations of 
O operator with respect to t such as Eq. (|23II26I) , we will 
derive another set of differential equations with respect 
to the initial value s, and it will be showed later that 
these equations are linear. Taking the derivative with 
respect to s of Eq. (IA10|) . one gets, 

^jz\{w\uMm 

= J2{z\{wm-^U-'[a,,H^,ilis)]Usm 

i 

= -iY,^{z\{w\Uta^{s)\0) 

i 

-i^\{z\{w\Ut[a,+i{s) + a,_i(s)]|0) 

i 

-iN{z\{w\UtY.f*e-^^^^d,{sm 

i 

- iN{z\{w\Ut ^ /i*e— '^e|(s)|0), (A12) 

i 

and noticing Eq. (|A4IIA7p . 

^jz\{w\uMsm 

= ~iuj{z\{w\UtL{s)\0) - 2iX{z\{w\UtL{s)\0) 

[ ds'ai{s,s'){z\{w\UtL{s')\0) + Nwl{z\{w\Ut\0) 

(A13) 



N / ds'a2{s',s){z\{w\UtL{s')\0). 



So, 



Oi = -iujOi - 2iXOi 

-N [ ds'ai{s,s')Oi{t,s',z*,w*) 



ds 



+ Nwl-N J ds'a2{s',s)Oi{t,s',z*,w*). 

(A14) 

Similarly, we can obtain 
d 



ds 



O2 = ioj02 + 2i\02 



+ N ds'a2{s,s')02{t,s\z* ,w*) 



Nz* + N / ds'ai{s' ,s)02{t,s' ,z* ,w*). 



Thus, we have obtained the differential equations for Oi 
and O2 with respect to s. 

Recall the master equation derived in Ref. [s^ . 



d_ 

di 



p = -i[Hs,p\ + [L, M{PtO\}] - [L\M{OiPt} 



[L\M{PtO\}]-[L,M{02Pt}] 



(A16) 



Define 



i?i(t, s) = Af{OiFj, R\{t, s) - M{PtO\}, (A17) 
R2{t, s) = M{02Pt}. Rlit, s) = M{PtOl}. (A18) 

then, from Eq. (jA14IIA15)) . we will obtain, 

d /"^ 
—Ri = -iujRi -~ 2i\Ri - N ds'aiis, s')Ri{t, s') 
"5 Jo 

+ N f ds'a*{s,s')Rl{t,s') 

-N I ds'a2{s' ,s)Ri{t,s'), (A19) 

J s 

with the initial condition 

Ri{t,t)^Lp. (A20) 
In the derivation above, the Novikov theorem 



(A21) 



M{wlPt}^ / ds'aiis, s')M{PtO\{t,s\z\w*)} 
Jo 

ds'aiis, s')R\{t,s'), 

10 

is used. 

Similarly, equation for R2 can be derived as 

— i?2 = iijjR2 + 2iAi?2 + N / ds'a2is, s')R2it, s') 
OS Jo 



N [ ds'aiis, s')R\it,s') 
Jo 



+ N ds'aiis', s)R2it,s'), 



(A15) 



(A22) 

with the initial condition 

R2it,t)=L^p. (A23) 
The solutions for R2 should take the following forms 
i?i = ci(t, s) ^ aip + C2 (t,s)^pa„ (A24) 

i i 

R2 = csit, s) V alp + C4(i, s) V pal (A25) 
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Substituting Eq. (jA24IIA25p into Eq. (|A19IIA22p . the 

equations for the cocfBcicnts can be derived as 

—ci = —iuJCi — 2iXci—N / ds' ai{s, s')ci{t, s') 
OS Jq 



+ N ds'al{s,s')c4{t,s') 



-N ds'a2is',s)ci{t,s'), 



(A26) 



—C2 = —iuJC2 — 2iXc2 — N / ds' ai{s, s')c2{t, s') 

OS Jq 



+ N ds'a2{s,s')c3{t,s') 



N / ds'a2{s',s)c2{t,s'), 



— C4 = iiujC4 + 2zAc4 + iV / ds' a2{sT s')c4^{t, s') 

OS Jq 



N / ds'al{s,s')ci{t,s') 



+ N ds'ai{s',s)c4{t,s'), 



with the initial conditions 





Cl{t,t) = 


= 1 




C2it,t) = 


= 


(A27) 


C3{t,t) -- 


= 1 




Ci{t,t) = 


= 



(A29) 



(A30) 



— C3 = iwc3 + 2iAc3 + iV / ds' a2{s,s')c3{t,s') 

OS Jq 



N / ds'al{s,s')c2{t,s') 



+ N ds'ai{s',s)c3{t,s'), 



(A28) 



Then, substituting Eq. (|A24IIA25|) into Eq. (UTH)) . the 
final master equation will be derived as it is shown in 
Eq. ([57]) . Note that the coefficients Ci{t) are determined 
by Ci{t,s) {i = 1,2,3,4) and the correlation functions 
aj(t, s) {j — I, 2) as shown in Eq. 
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